function res = trapezoidal_onestep(u,N,x_max,k)
Q = cal_Q(N,x_max);
I = eye(N);
res = (I-Q*(2/k))\(I+Q*(2/k))*u;